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O Abstract 

■ We study the properties of marginal distributions-projections of the phase 

space representation of a physical system-under relativistic transforms. We con- 
sider the Galileo case as well as the Lorentz transforms exploiting the relativistic 
^ • oscillator model used for describing the mass spectrum of elementary particles. 
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: 1 Introduction 



The concept of phase space arises naturally from the Hamiltonian formulation of clas- 
sical mechanics, and there have always been considerable efforts to give phase space 
picture of quantum mechanics too. Much of the thrust of these attempts lies in their 
ability to exploit classical analogues. Using these techniques, such as P-representation 
of Glauber and Sudarshan , the Wigner representation and the Husimi represen- 
tation [^, some quantum systems can be reduced to non-operator systems. However, 
the essential quantum nature of the problem is present in terms of the interpretation 
of the (apparently) classical variables. 
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Moreover, the so called quasi-probability distributions do not have the character- 
istics of classical probability distributions. Instead, by projecting the quasi-probability 
in a certain phase subspace, it is possible to obtain a genuine probability. In particular 
the projection of the Wigner function onto a straight line of the phase space was called 
'marginalization' procedure, and the obtained distribution 'marginal' 0, ^. 

Recentely, there has been a renewed interest on these marginal probabilities in con- 
nection with the tomographic imaging of a quantum state . Along this approach the 
marginals represent the shadows from which the state (or its phase space representa- 
tion) is reconstructed |^. 

The aim of the present paper is to study the properties of these marginals under 
relativistic transforms. The case of Galileo transforms results almost trivial, while for 
the Lorentz one we needed of a model having a covariant phase space picture. To 
this end we have studied the relativistic oscillator model used for describing the mass 
spectrum of elementary particles 

2 The "marginal" distributions 

Referring to the standard definitions given in the literature P, |^, by 'marginalization' 
one should mean a line integral in the phase spce {q,p} of the Wigner function W{q,p), 
i.e. 



where 6 is the angle orientation of the line, w becomes a probability distribution for 
the variable x, depending parametrically on 6. 

One can go beyond this definition ; let us consider the phase space transformation 
as a generic linear combination of position q and momentum p 




(1) 



X 



fiq + vp 



(2) 



p 



V 



/i'g + v'p 



(3) 
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and consider it as a real symplectic transformation belonging to the group Sp{2,R), 
i.e. 

A.A^^.. A^(;, :,), .^(_«^ J). (4) 

Then, we can consider the 'projection' along the phase subspace characterized by 
the transformation (|^) 

w{X;n,i/) = J W{q,p)6{X — nq — vp) dq dp , (5) 

which can be intended as a marginal distribution too. Of course, Eq.(|l|) represents a 
particular case of Eq.(^ whenever a mere rotation in the phase space is considered. 

The above definition could also be extended to phase spaces of higher dimensions. 
For example, in the case of two-dimensional system we will have a phase space {q = 
ill: 12), P = {pi,P2)}, hence we can intruduce the transform 

(X,P) = A(g,p)^ (6) 



where now A is a 4 x 4 real symplectic matrix, i.e. 

1 0\ 
n 1 

(7) 



AaA^ = cr; a 



f 1 ^^ 
00 01 
-10 00 
V -10 0/ 



The components Xi, X2, Vi,V2 are related to the homogeneos symplectic group 5*^(4, R). 
In particular, X has the following components 

Xi = flq + up, (8) 
X2 =pi'q + v'p, (9) 

with /I = (An,Ai2); z7 = (Ai3,Ai4); pi' = (A2i,A22); z^' = (A23,A24). 

Therefore, a marginal distribution w{X; fl, z7, /I', z7'), may be introduced as a prob- 
ability distribution for the variable X, with a dependence upon the parameters char- 
acterizing the matrix A 

w{X; jl, V, jl', u) = J dq dpW{q, p) 5{Xi — jlq — up) 



X 5{X2-p^q-Vp). (10) 

3 Marginal distributions and relativity 

Now, in order to study the properties of marginal distributions under relativistic trans- 
formations we need to know the transformation properties of the wave function of the 
system. 

At first it is instructive to consider the Galilei transforms 

q = q — vt] p' = p — v; t' = t , (11) 

where we have considered a particle of unit mass and v represents the relative velocity 
between the two reference frames. 

Since in this case the wave function transforms, independently of the assumed 
model, as 

^(g, t) ^ ^'{q, t) = exp (^tvq - ijt^ ^{q - vt, t) , (12) 
then, the Wigner function correspondingly transforms as 

W{q, p, t) W{q\ p\ t) = W{q -vt,p- v, t) . (13) 

As consequence of (0) we immediatly obtain 

w{x, /i, u) — i> Wy{x, /i, z/) = Wo{x — fiVt — UV, fi, u) , (14) 

from which a simple shift in the distribution, in moving from one reference frame to 
the other, results. 

Instead, in the case of Lorentz transformations we have to consider a specific model. 
3.1 The relativistic harmonic oscillator 

Let us consider the relativistic oscillator model introduced by Feynmann et al. to 
describe a hadron consisting of two quarks bound togheter by a harmonic oscillator 
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potential of unit strength 




where Xa and Xf, are space-time coordinates for the first and second quark respectivly 
(we are using natural units, h = c = 1). This partial differential equation has many 
different solutions depending on the choice of variables and boundary conditions. Here 
we follow the treatment of Ref . [0 . 

In order to simplify the Eq. (p!5|), let us introduce new coordinate variables 



The four- vector X specifies where the hadron is located in space-time, while the variable 
X measures the space-time separation between the quarks. In terms of these variables 
Eq. (|15D can be written as 



X = {Xa + Xb)/2- 



X = {Xa- Xb)/2 . 



(16) 




(17) 



This equation is separable in the X and x variables. Thus 



(f){Xa,Xb) = f{X)%l){x) , 



(18) 



and f{X) and il){x) satisfy the following differential equations respectively 




(19) 



(20) 



Equation (|T9|) is a Klein-Gordon equation, and its solution takes the form 



/(X)=exp[±zP^X^] 



(21) 



with 



-P2 = - 



P^P^ = =ml + {X + l) 



(22) 
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where M and P are the mass and four-momentum of the hadron respectively. The 
eigenvalue A is determined from the solution of Eq. (^OD . 

As for the four momenta of the quarks Pa and pf,, we can combine them into the 
total four-momentum and momentum-energy separation between the quarks 

P = Pa+Pi>] p=V2{pa-pb). (23) 

P is the hadronic four-momentum conjugate to X. The internal momentum-energy 
separation p is conjugate to x provided that there exist wave functions which can be 
Fourier transformed. 



The four-dimensional equation (|20|) is separable in at least thirty-four different co- 



ordinate systems Since we are quite familiar with the three-dimensional harmonic 
oscillator equation from nonrelativistic quantum mechanics, we are naturally led to 
consider the separation of the space and time variables, and write the equation ( pOD as 

+ ^ + - ^']) ^(^) = + l)^(^) • (24) 

If the hadron moves along the Z direction which is also the z direction, then the 
hadronic factor /(X) is Lorentz-transformed in the same manner of a scalar field. The 
Lorentz transformation of the internal coordinates from the laboratory frame to the 
hadronic rest frame takes the form 

x' = X, y' = y, 

z' = iz-pt)/{l-PY^\ (25) 

where /3 is the velocity of the hadron moving along the z direction. The primed 
quantities are the coordinate variables in the hadronic rest frame. In terms of the 
primed variables the oscillator differential equation is 

+ ^ + [x'^ - t"]\ ^{x) = (A + mx) . (26) 



This form is identical to that of Eq. (0), due to the fact that the oscillator differential 
equation is Lorentz-invariant . 

Among many possible solutions of the above differential equation, let us consider 
the form 

1 \ / 1 \ (''+^+"+^)/2 / I \ 1/2 



X exp 



a\b\n\k\ 



Ha{x')H,{y')Hn{z')Hk{t') 



--[x +y +z +t ) 



(27) 



where a, b, n and k are integers, and Ha{x'), Hh{y') . . . are the Hermite polynomials. 
This wave function is normalizable, but the eigenvalue takes the values 



\ = {a + b + n — k). 



(2^ 



Thus for a given value of A, there are infinitely many possible combinations of a, b, 
n and k. The most general solution of the oscillator differential equation is infinitely 
degenerate []r3|. The simplest way to avoid this problem (at least to render finite the 
degeneracy), is to invoke the restriction that there should not be time-like oscillations 
in the Lorentz frame in which the hadron is at rest, and that the integer k in Eqs. ( pTf) 
and (^Sf) be zero |13|, In doing so we are led to the question of maintaining the 
Lorentz covariance with this condition. 

When the hadron moves along the z axis, the k = condition is equivalent to 



The most general form of the above condition is 

d 



Pf,[x^ + 



dx, 



l/jf}{x) 



0. 



(29) 



(30) 



Thus the k = condition is covariant. Once this condition is set, we can write the 
wave function belonging to this finite set as 

■ 1 \ / 1 \ ("+^+")/2 / 1 \ 1/2 



TT 



a\b\n\ 



H,{x')H,{y')H^{z') 



X exp 



/ /2 . /2 I /2 I i/2\ 

--{x +y +z +t ) 



(31) 



Except for the Gaussian factor in the t' variable, the above expression is the wave 
function for the three-dimensional isotropic harmonic oscillator. 

Since the above oscillator wave functions are separable in the Cartesian coordinate 
system, and since the transverse coordinate variables are not affected by the boost 
along the z direction, we can omit the factors depending on the x and y variables when 
studying their Lorentz transformation properties. Hence, the solutions satisfying the 
subsidary condition (|30D take the simple form 



/ 1 X 1/2 



1 



2 , ,/2x 



(32) 



with \ = n. This normalizable wave function, without excitations along the t' axis, 
describes the internal space-time structure of the hadron moving along the z direction 
with the velocity parameter (3. If /5 = 0, then the wave function becomes 
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(33) 



Thus 

r^{z,t) = r,{z\t'). (34) 

We have therefore obtained the Lorentz-boosted wave function by making a passive 
coordinate transformation on the z and t coordinate variables. 

3.2 Covariant phase space 

It is possible to construct a covariant phase space for the relativistic harmonic oscillator 
by following Ref. [l^. Let us consider at first the Gaussian factor of the wave function 



(|33|), which practically corresponds to the ground state, 

^°(^,t) = (^)"\^v{-{z'+t')/2) , (35) 
and introduce the light cone coordinates 

u = {z + t)/V2, v = {z-t)/V2. (36) 



The latter transform as 

It is easy to see that the product uv is Lorentz invariant. By using such coordinates, 
the wave function ( |55D can be rewritten as 

^^,{z,t) = i;',iu,v) = (^)'^'exp {^-{u^ + v^) /2) , (38) 

and, if the system is boosted, it becomes 



exp|-^-J^^«^ + ^„^J|. (39) 

Practically, it undergoes a continous deformation as f3 increases. 

Analogously to the Eq. (|3^), we may define for the momentum and energy the 
variables 

Pu = iPz-Po)/\^, Pv = iPz+Po)/\^, (40) 
and the momentum-energy wave function will be given by 

0S(P.,P.) = (^)/^S(;^,t)e-^(^^-*^°)dzdt, (41) 



Hence, we deal a four dimensional phase space where the Wigner 

function can be defined in a canonical way 

W^{u,pu]v,p^) = (^^^ J (V'°(M + x,f + ?/))* ?/^J(M-x,t;- y) 

X exp [2i {puX + PvV)] dxdy . (43) 

After the evaluation of the integral we obtain 



l^ ('i+(3.2 1-." 2 



which is manifestly co variant. 



3.3 The properties of marginal probabilities 

Having a four dimensional (covariant) phase space, the marginal distributions can be 
defined analogously to the case of Eg. (p^) . Practically, a marginal distribution will be 
a projection on the plane {W, V} determined by the equations 

U = fiiU + fi2Pv + ^iV + l'2Pu , (45) 

V = Ciu + C2PV + Viv + V2Pu , (46) 

where we have /Xj, z/j, Q, rji eH, (z = 1, 2). 

Without lost of generality we do not specify the constraints on these parameters, 
since they will be related to the space-time asimmetry in the commutation relations, 
which is an hard problem to face in making the relativistic quantum mechanics, and 
goes beyond the scope of the present paper. 

Then, we define 



Wf3{h{,V;a) = J dudvdpudpyW^{u,pu;v,pv) 

X 6{l{ - jJ,iU - ^2Pv - T^lV - l'2Pu) 

X 5 (V - CiM - C2PV - Viv - mPu) , (47) 

where a = {/ij, z/j, Q, rji}, {i = 1,2). As limiting cases we have /ii = r^i = 1, and all the 
other parameters equal to zero, then 

wp{U,V;a) = \ij^^{u,v)\' ; (48) 

or otherwise, for 1^2 = C2 = ^, and all the other parameters zero, then 



wp{U,V;a)= <Pl{pu,Pv)' ■ (49) 



2 



It is clear from the Eq. (|47|) that the marginal distribution, with such definition, is 
not covariant, but we may rescale the variables, due to the Wigner function covariance. 
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to have 



wJ(W,V;(t) = j dudv'dp'^dp^Wo{u',p'^;v',p'^) 

(5 (V - CiM' - VlV' - V2P'u - C2P'v) 



X 



where we take into account the invariance of the measure and we set 



(3 



1 + P\ 

From Eq. (|50| ) it immediately follows 



1 + P 



1/2 



1/2 



.1-/3 
1-P 



nl/2 



(50) 



(51) 



1/2 



1,2) 



ii;°(W,V;a)=u;S(W,V;a^), (52) 

where indicates the parameters (|5TD. 

Eq. ( p2D defines the transformation properties of the marginal distributions; prac- 
tically we get that different marginals correspond to the same measurement, but per- 
formed in different frames. This means that the boosts connecting several frames could 
be useful to vary the parameters characterizing the marginal distribution. This inter- 



pretation is in agreement with that given in Ref. ||T6[ for the two-mode 'symplectic 
tomography' . 

The above results can be extended to the excited states of the relativistic oscil- 
lator as well; the only additive factor one has to consider is the Hermite polynomial 
multypling the Gaussian of the ground state. 

Of course, the result of Galilei transforms cannot be obtained as a limiting case of 
the Lorentz transforms. 



4 Conclusions 



In conclusion we have studied the properties of the marginal distributions under rel- 
ativistic transformations. Since they contain all the information about the quantum 
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state of a system, other probabilities related to different observables, could be derived 
from them, as well as their properties. 

Finally, the discussed properties do not concern only fundamental questions, but 
could become interesting in quantum optics where "optical mesons" enter in the reality 



T7| , and in particle physics where, hopefully, the quantum state tomography concept 
could be applied. Since by repeated measurements one can build up the marginal 
probabilities, one is lead to ask the following question: which observables should be 
measured to this goal for example in high energy processes? This subject will be 
addressed in future works. 
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